> 

On 



email: bergere@spht.saclay.cea.fr 



1 Introduction 



THE PARTITION FUNCTION FOR "COMPOSITE PARTICLES" 

by 

M. C. Bergere 

^ - Service de Physique Theorique, CEA-Saclay 

Q\ ■ F-91191 Gif sur Yvette cedex, France 

Oh! 

in 

Abstract: We calculate the partition function for "composite particles". 
For any finite number of states d, and in the following two cases: 1°) all states 
have the same energy, 2°) the energy is linearly distributed over the states, we 
transform the partition function into a finite sum of terms which exhibit trivially 
their d dependance. The infinite d thermodynamic limit is obtained as well as 
the finite <i-size corrections, etc... In the second case, we obtain the finite <i-size 
corrections to "the universal chiral partition function for exclusion statistics". 

G\ 

a 

Fractional statistic was introduced in the literature in order to explain various 
effects encountered in the physics of condensed matter such as the fractional 
quantum Hall effect, the high temperature superconductivity, the properties of 
the spin 1/2 antiferromagnetic chains... and was explicitely observed in several 
theoretical models like for instance the two-dimensional anyon model' 1 ' or the one- 
dimensional Calogero-Sutherland model t 2 '. This statistic generalizes the exclusion 
Pauli principle and interpolates between bosons and fermions statistics in the 
sense that it takes g states to add one more particle to the system (g = for 
bosons and g = 1 for fermions) where g is fractional (g = l/m); Haldane^ 
proposed that in the fractional statistic, the number of possible configurations 
for N particles on d states should be 

Of course, this formula makes sense only statistically for large N and d, otherwise 
[d + (1 — g) (N — 1)] is a fraction and the combinatoric C symbol calculated with 
the r Euler function gives a fractional number of configurations. However, we 
note that for finite N of the form N = pm + 1, the combinatoric C symbol 
makes sense and gives an integer number of configurations. This remark leads 
us to propose^ a specific definition of the "composite particles" which naturally 
generates the fractional statistic: 
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A "composite particle" is a set of m particles and d > I states with the 
constraint that the (I — 1) bottom states are empty and the I state (from the 
bottom) has at least one particle. We define an uncomplete " composite particle" 
as a set of < r < m particles with no constraints on the empty states. 

The main constraint in this construction is that a non empty state is entirely 
included inside a " composite particle" and cannot be splitted over several " com- 
posite particles" (fig.l). The number of possible configurations for the above 
definition of "composite particles" has been found in W as 

C d+N-l-lE(^) ( 2 ) 

where E (x) is the integer part of x. If N = pm + 1, this formula gives the same 
result as Haldane's formula, but (2) makes sense for any N = pm+r (0 < r < m). 
For large N, the value of r is irrelevant and the asymptotic properties of (1) 
and (2) are equivalent. In that sense, the "composite particles" reproduce the 
fractional statistic introduced by Haldane and extend it naturally to any finite 
iV and d; in addition, it gives a geometrical interpretation of this statistic. 

If all states have the same energy E, it is found that for large iV and d 
( lg (N\) ~ iVlg (N/e) ), the average number of particles per state is given by ^ 

n (y, g) = w t 1 x . < - (3) 

W{y,g) + g g 
where y = exp (—E/kT) and where the function W (y,g) satisfies 

[W(y,g)] g [l + W(y,g)t B =y' 1 (4) 
In appendix A, we show that 

w («, 9 ) = g - + 9 - 1 + (, - 1) t - '" 9 " t+iY. ~ n + l) { - y) " (5) 



n=l 

as already found by Ramanujan' 6 '. 



More recently' 7 ', it has been proved that the "composite particles" with iV 
particles were in a one-to-one correspondance with the Young tableaux with at 
most iV rows. As a consequence, and using the Bethe Ansatz, it was possible to 
reproduce from the "composite particles", the eigenstates of the Hamiltonians for 
the models of Calogero-Sutherland' 2 ] and of Ruijsenaars-Schneidert 8 '. From the 
complete knowledge of the spectral properties of these models, it was then possi- 
ble to define a momentum space representation for the "composite particles" in 
terms of creation operators attached to some Young tableaux; these creation op- 
erators are generated from the normal ordered product of vertex operators when 
expanded in terms of Jack or Macdonald' 9 ' homogeneous, symmetric polynomials. 
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Similarly, a position space representation has been proposed for the "composite 
particles" in terms of appropriate products of vertex operators. 

In this publication, we calculate the partition function for the "composite 
particles" and for any finite number of states d, in the following two cases: 
1°) when all the states have the same energy E, we find 

Sup(m,l)— 1 

Z(y,l,m)= a j(vJ^ m ) [/3j(y,l,m)] d (6) 

where y = exp (—E/kT) . We determine the <i-independant functions Oij (y,l,m) 
and j3j (y, I, m) ; this result generalizes the bosons partition function (1 — y)~ d , 
and the fermions partition function (1 + y) d . 

2°) when the energy is linearly distributed over the states, that is when the 
state i contributes to the partition function by a factor 

x^y = exp [- (E+(i- 1) v) /kT] v>0 (7) 

In that case, we find 



Sup(l,m) — 1 

Z(x,y,l,m) = Y a j(yi x i l i m ) 



d-l 



Yl^j (x l y,x,l,m) 



i=0 



ij {x d y, l, m) 

(8) 



and we determine the d independant functions aj(y,x,l,m), (3j (y,x,l,m) and 
7j (y, x, I, m) . The case 2°) is essentially a a;-deformation of the case 1°). In both 
cases, we transform the partition function for any finite d, into a finite sum 
of terms which exhibits trivially their d dependance; consequently, our result 
provides naturally the large d asymptotic behaviour, the finite d— size corrections, 
etc... of the partition function. 

In the litterature, the function (3 Q (y,l,m) in (6) which describes the large d 
thermodynamic limit is known^ as 

/3 (y,l,m) = l + V(y,l/m) (9) 

where V (y,l/m) = [W (y, g)] 1 satisfies the equation 

v{i + vy- x = y (10) 

On the other hand, we define the thermodynamic limit in case 2°) as the limit 
x = £rf _> i anc i d _ > we obtain for the leading contribution 

d—l f A \ 

hW^oo J[P (z*y,Zv,l,rn) ~ e *p (^-^ [B (y) - B (11) 



3 



where 



B(y) = [ y ^]n[l + V(v,g)} 
Jo v 

V(y,g) 



B{y) 



l + V(y,g) 



+ Uny. In [1 + V (y,g)] (12) 



and where L (y) is the Rogers dilogarithm function. The Rogers dilogarithm 
function and dilogarithm identities were previously used by Nahm, Recknagel and 
Terhoeven' 10 ' to calculate the structure constants and, in particular, the effective 
central charge of certain conformal field theories. Later, Hikami' 11 ' calculated the 
effective central charge for an ideal g-on gas as 

-S/f-^-'M-lf (13) 

where £ H satisfies the equation 

?H = 1 " ^ ( 14 ) 

More recently, the function B (y) was obtained by Berkovich and McCoy' 12 ' us- 
ing the steepest descent method in order to evaluate the thermodynamic limit 
of what they call "the universal chiral partition function for exclusion statis- 
tics". Their partition function was defined in the context of conformal field theory 
and of Rogers- Ramanuj an identities. They demonstrate the equality between 
their partition function and Haldane's exclusion statistic. "The universal chiral 
partition function for exclusion statistics" has been defined for several different 
kinds of particles where g is replaced by a matrix Gab] in Ref [13] several authors 
have explored the properties of these statistics in order to apply them to physical 
situations in condensed matter physics. 

In section 2, we calculate the partition function when all the states have the 
same energy. In section 3, we give the general expression of the partition function, 
when the state i contributes to a factor Xi = exp (—Ei/kT) , as a sum of non- 
symmetric polynomials of the x^s. In section 4, we apply section 3 to the case 
where the energy is linearly distributed over the states as in equation (7); in that 
case, the partition function can be written as the following sum 

Z(x,y,l,m) = l+ E xi ^ m+Pr)l C Z\™-i) + r-i (*) (15) 

p=0 r=l 

It must be noted at that point that "the universal chiral partition function for 
exclusion statistics" as considered by Berkovich and McCoy' 12 ] and which orig- 
inated in the work of Kedem, Klassen, Mc Coy and Melzer' 12 ] in the context 



4 



of conformal field theory, does contain the symbol C^ 1 _ g ^ N _ 1 ~ ) (x) (generalized 
to several kind of particles with a matrix g); however, their calculation of the 
thermodynamic limit is made for one kind of particles at infinite d so that the 
symbol C^ (1 _ 9)(iV _ 1) (x) is replaced by [(x, x)^ 1 where (x,x) n = fliLi i 1 ~ x *) 
and their limit when x — > 1 is a special case of ours. 

In the rest of section 4, we transform (15) into (8) for any finite number d of 
states. 

Let us conclude this introduction by mentionning the fact that our construc- 
tion of the "composite particles", for any finite number of particles N and for 
any finite number of states d, reproduces perfectly well previous results known in 
the literature, provides finite ci-size corrections to these results and, maybe more 
important, gives a quasi geometrical interpretation to fractional statistic. 



2 The partition function Z(y, /, m) 

The partition function 

Z(y,l,m)= J2 Ci,^) y N ( 16 ) 

N=0 

may also be written in terms of the number of " composite particles" p and of the 
remaining number of particles r in the uncomplete " composite particle" : 

*(<¥) m 

Z(y,l,m) = l+ £ £ y^ (17) 

p=0 r=l 

The sums in (17) can be calculated for any finite number of states d. The purpose 
of this section is to transform the polynomial in y which defines the partition 
function into a finite expansion of the type 

Sup(l,m)— 1 

z(y,i,m)= a j(yJ' m ) [Pj(yJ, m )] d (18) 

3=0 

where the sum over j is finite and where the functions aj (y, I, m) and /3j (y, I, m) 
are d independant. The easiest examples are given by the cases of the bosons 
and the fermions statistics where 



oo 1 

Zb = Y1 O-i y N = Tr-- d (19) 
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d 

z F = J2 c?y N = {i + y) d (20) 

N=0 

In the general case, the structure of the function Z(y,l,m) is twofold. A first 
effect can be traced in the special case m = 1 and results in the introduction of 
functions related to W(y,g) defined in (5). A second effect appears when m > 1 
and is expressable as a sum of similar expressions which differ by the choice of a 
m th root of unity. This second effect becomes particularly simple for I = m since 
in that case the function W(y,g = 1) = y -1 . 

Let us first consider the second effect that is the case of "composite fermions" 
where / = m > or g — 1. We define 

Z r (y,m,m)=5 r , m + £ C^ r + _\ y? m+r (21) 
An easy calculation gives 

1 m— 1 

Z r (y, m, m ) = -J2 9"™ (1 + g"!/)**" 1 (22) 

n=0 

where g = exp {2m /m) . Finally, summing over r from 1 to m, we get the partition 
function for the "composite fermions" as 

Z (y, m,m) = - Y (1 + ^ ~ \ (1 + q n yf (23) 

which is of the form (18). 

The thermodynamic limit, at large d is dominated by the root of unity equal 
to 1 that is n = 0; we get 

Z th (y, m, m) = - Q + vFjlL (1 + y f (24) 
m y 

which shows that we recover the usual fermion partition function corrected by m 
dependant "finite d— size effects". At zero energy, this correcting factor is simply 
(2 m - l)/m. 

The first effect (for m — 1) is obtained when calculating 
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This sum can be calculated either from a recurrence on Z as a function of d (like 
in section 4.3 and appendix B) or by introducing a contour integral representation 
of the combinatoric symbols as already presented in Ref. [14] 

Ci = ^-idz^f (26) 

where the contour Co surrounds the point z — 0. It is useful to note that if the 
contour surrounds z = and z — — 1, then the contour integral gives zero for 
N > E so that the geometrical series in N can be summed to infinity 

provided that the series converge on the contour of integration. We obtain for 
I > 

ZiyJA) = 2^f dZ hlV />0 (27) 
mk j Cl z(l + z) -y 

The convergence of the series requires \z\ |1 + z^ 1 > y which is the case for large 
\z\ if I > (we remind that < y < 1). Consequently the contour C\ surrounds 
the / roots Vj (y, I) of the polynomial equation 

z{l + zf- 1 = y (28) 
If we expand the denominator in (27) around each pole z — Vj (y, I), we obtain 

which is again of the form (18). 

Several remarks should be made at this stage: 

1°) the I roots Vj (y, I) are independant of the number of states d, 

2°) when y — > 0, the root V (y,l) ~ y; the (I — 1) other roots — > — 1 + 

O (y~^j ■ All the roots can be expressed in terms of a unique function V (y, I) 
defined in appendix A: 



V (y,l) = V(y,l) 

Vj(y,l) = -l-y (V+'^T^l) J = 1,-,J-1 (30) 

where Q = exp (in/ (I — 1)) , 

3°) the sum over all roots is such that Z (y, 1, 1) is a polynomial in y of degree 
E { d+ \~ 1 ) ; the roots Vj (y, I) for j > contribute to powers of y > E ( ^j" 1 ) , 

4°) from equation (28) and the zero y limit, we see that V (y, I) = [W (y, Z)] -1 
defined in (5). 
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In the thermodynamic large d limit, the polynomial in y which defines the 
partition function in (25), becomes an infinite series which can be summed into 

[l + V(y,l)] d+l 

Ztt(y,U)= 1 + lvM l >° (3 1 ) 

expressed in terms of a unique function V (y, a) defined in appendix A:This result 
is also valid for I — where V (y, I) — yj (1 — y) . 

We now write the expression for the partition function Z (y, I, m) defined 
by the sum (17) for any I and m. This partition function contains both effects 
described above. The details of the proof are given in appendix A; we get 

Sup(m,l)-1 , w , T/\d+- 

Z(y,l,m)= E + i+V i / m (32) 

where q = exp {2m /m) and where k (j) is an integer equal to j for j = 0, ...,m — 1 
but depends upon the definition of (1 + Vj) 1 ^ in (A17) for I > m and j = 
m, — 1. The functions Vj are the roots of the polynomial equations 

z m - (1 + z) m ~ l y m = m>l (33) 

z m (1 + z) l ~ m - y m = I > m (34) 

and may be written in terms of a unique function V (y, a) defined in appendix 
A. 

In the thermodynamic large d limit, the partition function is dominated (at 
least for y small enough) by the root V 0) so that 

Zth iy , m) = (1 + VO./M));-! (35) 

(l + V(y,l/m))™ - 1 rn + lV{y,l/m) 
The average number of particles is easily obtained from 

N^y% (36) 

and 

dV(y,l/m) = V(y,l/m)(l + V(y,l/m)) 
V dy 1 + ± V (y,l/m) 1 } 

In the large d limit, we obtain the average number of particles per state as 

n(y,g) = hm d^ oo — = — — 38 

d l + gV(y,g) 

which shows again from (3) that V (y, g) = [W (y, g)] 1 . 
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3 The general structure of the partition func- 
tion 



We now construct the partition function when the state i has a given energy E t 
for i — 1, d. We introduce the reduced variables Xj = exp (—Ei/kT) . We first 
consider the uncomplete "composite particle" with r particles on S states. Its 
contribution to the partition function is 

y r x i- x 8 S (39) 

{^>o} 

where rj is the occupation number of the state % and y is a variable which counts 
the number of particles. Of course the sum in (39) is restricted to Ylt=i r 'i = r - 
This contribution can be written in terms of the symmetric, homogeneous (of 
degree r) Schur polynomial attached to the Young tableau with one row and r 
squares 



y r S r (x 1 ,...,x s ) (40) 

Now, a "composite particle" with S states has m particles on 5 — I + 1 states; 
moreover, we have to make sure that the state "<5 — I + 1" has at least one 
particle; consequently, the contribution of such a "composite particle" to the 
partition function is 

y m [S m (xi, xs-i+i) - S m (xi, x S -i)] = y m xs-i+i S m _i (x 1: ...,x S -i+i) (41) 

We now organize the p "composite particles" each of which having Si states 
(i = l,...,p) and the uncomplete "composite particle" with 5 p+ ± states. Let us 
introduce the notation 

k 

A k = J2^n (42) 
n=l 

Then, the partition function is 

z = i+ £ E y pm+r E(n(^ fc -m^-i(^_ 1+ i,...,x Afc _ m )) . 

p=0 r=l {Si} k=l 

■ S r (x Ap+1 ,...,x d )) (43) 

where the sum over {5i} is restricted by A p+ i = d where d is the total number 
of states. The sum over {Si} may be performed to obtain 

Z = l+ J2 ^Z N ( Xl ,...,x d ) y pm+r (44) 

p=0 r=l 
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where Zn (xi, ...,Xd) is a homogeneous (of degree N = pm + r), non symmetric 
polynomial of the variables Xi : 

z N ( Xl ,...,x d )= Yl x i- x d ( 45 ) 

{pi}eA N 

,m,( 

In (45), the set Ajv m / of possible p^s is defined by three constraints: 
1°) 

d 

J2Pi = N (46) 



i=i 



2°) let 



7T i = ^2p j i = l,...,d (47) 



then, 



{m,2m,...,pm} C {m} (48) 
3°) we consider the indices (ii, ...,i p ) such that n iq = qm, then 

Pi q +i =Pi g +2 = ••• =Pi q +i-i = g = l,...,p (49) 
The number of monomials in Z N is / w i \ ■ 

d+N-l-lE{^=±) 

4 The partition function when the energy is lin- 
early distributed over the states. 



In this section, we apply the structure obtained in section 3 to the case where 
the energy for the state i is 

Ei = E+(i-l)v v>0, i = l,...,d (50) 

The reduced variables are defined as x l ~ l y = exp (—Ei/kT) . The Schur polyno- 
mials in (43) are known to be 



S n (l, x, x 2 , x x ) — C£ +n _ 1 (x) (51) 
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where 

CI (x) = [ X \ X) \ (52) 

and 

(a; x) n = (1 - a) (1 - ax) ... (l - ax^ 1 ) (53) 

Using the homogeneity property of the Schur polynomials, we obtain for the 
partition function 

E{— r~ ) m p 

Z(x,y,l,m) = l+ J2 E^E^II^^-i^) • 

p=0 r=l {5,} i=l 

(54) 

where 

P-i 

i/ = m ^ Ai + (r + 1) A p - pi (55) 
i=i 

and Aj is given in (42). Next, we must sum over all number of states (Si > I if 
1 < i < p, and S p+ i > 1) with the restriction that A p+1 = d. To perform this sum, 
we apply the usual technique of the generating functional; we multiply Y2{Si} ■■■ 
in (54) by 

zt-K.J^zllT 1 (56) 

we sum all S^s from I to oo (i — 1, ...,p) and S p+ i from 1 to oo and we apply the 
relation 



5=0 

We obtain for J2{s t } •••^ n (54) 



E^W*)^^ (57) 



x (*?W)i TT 1 — (58) 

where = (p — 1) m + r + 1. If we set all z^s equal to z in (58), we get 

x (*fW)l 1 (5Q) 
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Finally, the partition function is obtained by determining the term of (59) in 
z d- P i-i _ Tj s i n g (57) backwards, we get 

*(<¥) rn 

Z(x,y,l,m) = l + £ E ^ m+pr)l CZ + ( L-i )+r -i (*) iT^ (60) 

p=0 r=l 

This expression is "the universal chiral partition function" for one kind of parti- 
cles. 

In the remaining part of section 4, we compute the sums over p and r in (60) 
for any x, y and finite d in the same way as in section 2. We first remind to the 
reader the special cases of the bosons and of the fermions; then, we distinguish, 
like in section 2, the case I — m (g — 1) and the case m — 1, before writing the 
general result which is proved in Appendix B. 



4.1 Bosons and Fermions 

The partition function for bosons is given by the expansion (60) with / = 

oo 

Zb = Y1 C&n-i (*) y N = jz77v (61) 
n=o \y^ x )d 

Similarly, the fermions partition function is given by (60) with I = m = 1 

z F = Y.x N -^c»{x)y N = (-v;x) d (62) 

N=0 

Clearly, for x — 1, we recover the well known results of section 2. The large 
d limit for x < 1 is given by [(y;^)^] 1 for bosons and (—y^x)^ for fermions. 
These finite results (which are not exponent ionnally large with d) imply that 
the average number of particles iV defined in (36) is a finite quantity, so that 
the average number of particles per state n(x, y, g) defined in (38) is zero. An 
interesting thermodynamic limit is obtained when x — > 1 and d — > oo. Let us 
write x = ^ so that (7) becomes 



i-l 



C d y = exp 



- ( E + % —± (£oo - E)^j /kT 



(63) 



At large d, the energy spectrum becomes linearly continuous from E to E^. Using 
the exponential representation 



1 
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= exp(-f^^f) (65) 

and performing a 1/d expansion, we obtain the thermodynamic limit for the 
bosons as 

= «•> E a - o $ + \ ± a - n $ + o (I) ) (66) 

\ ' n=l n=l / 



that is 



^B,th = J exp [L* 2 (y) - L* 2 (£y)]) + O Q) (67) 



where the dilogarithm function is 



Li 2 (y) = - / — ln(l-u) (68) 



Similarly, for the fermions 



\ ' n=l n=l / 



(69) 



that is 



where 



^ = exp [Li+ ( y) - Li+ &,)]) + O Q) (70) 

Li+ (y) = / - In (1 + u) = U 2 (y) - -Li 2 (y 2 ) (71) 
Jo u z 

We show in subsection 4.4 that the functions Li 2 (y) and Li + (y) can be expressed 
in terms of the Rogers dilogarithm function L(y). 

In the thermodynamic limit, the average number of particles per state is 

n B (e,y) = 1 ^=f (72) 



In 



1+2/ 



n F (Z,y) = ^%- (73) 

Of course, when £ — > 1, we recover the average number of particles per state of 
section 2. We remind that in formulas (67-70-72-73), y is related to E and £y 

2 2 

is related to E^. Finally, let us remind that Li 2 (1) = \cb and Li + (1) = \cf 
where cb = 1 is the central charge for bosons, and cp — 1/2 is the central charge 
for fermions. 
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4.2 "composite fermions" 



Now, we consider the case of "composite fermions" (I — m > or g — 1). We 
define 



P(P-1) — 2 



Z r (x, y, m, m) = 5 r , m + £ X ^^ +P rm {x) yPm+r (?4) 

p=0 

An easy calculation using (62) shows that 

r(m-r) m —l 

Z r (x, y,m,m) = X - V g' rw (-q^^y; x) (75) 



n=0 



where g = exp {2m /m) . The partition function Z (x,y,m,m) is obtained from 
(75) by summing r from 1 to m. 

The thermodynamic large ci-limit when x = £ d — > 1 can be performed simi- 
larly as for the fermion case. Again, the asymptotic behaviour of (75) is domi- 
nated by the term n = 0; using (53) and performing the ^ expansion, we obtain 

+0 (±) (76) 



We may sum the result over r from 1 to m: 



m 

Z tk (ty,m,m) = (j^P) '^mft, " (^|=T (</) ~ 

+0 (i) (77) 

This result shows that the leading exponential term is the same as for fermions 
(70), but here, we obtain some m dependant finite rf-size corrections. 



4.3 The partition function Z(x,y,l,l) 

The construction used in this subsection is also valid for x = 1 and provides 
another proof of the result of section 2 without using contour integrals. Below, 
we need to specify the number of states d. 



14 



The partition function is 



Z (x, y, 1, 1) = Z d (y, x)= C d + (i-D(N-i) (*) V* (78) 



Using the fact that the sum contributes nothing for N > E 1 ) because of 
the combinatoric symbol C, and the fact that 

CI (x) - Cl_, (x) = CtX (x) (79) 

we have 



Af(iV-l)i 



Z d (y, x) - Z d ^ (xy, x) = x 2 C d+(i-i)(N-i)-i ( x ) V* ( 80 ) 



N=l 



Consequently, shifting N by one unity in (80), we just proved the following re- 
cursion relation 

Z d (y, x) = Z d _i (xy, x)+y Z d ^i (x l y, x) (81) 

This equation (which is also valid for / = or 1) makes sense for any number of 
states d down to zero provided that we introduce (according to (78)) 

Z-i(y,x) = Z_ 2 (y,x) = ... = Z_ t+1 (y,x) = 1 

Z-i(y,x) = (82) 

Now, we suppose / > 1 and we define 

i-i 

V d (y, x) = Z d (y, x ) + J2 W i (y, x) Z d-i (x'y, x) (83) 

i=i 

The unknown d independant functions Wi (y, x) are fixed by the existence of a d 
independant function F (y, x) such that 

d-i 

^ d (y, x) = [l + F (y, x)} ^ (xy, x) = ]J[l + F (x'y, x)] ^ {x d y, x) (84) 

i=0 

The constraints on the functions Wi (y, x) which ensure (83) and (84) are 

W l (y,x) = F(y,x) (85) 
[l + F(y,x)]Wi(xy,x) = W i+1 (y,x) i = 1, Z — 2 (86) 
[1 + F(y,x)] Wi- X {xy,x) = y (87) 
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which means that 

fe-2 



W k (y,x) = Y[[l + F(x i y,x)] F (x fc_1 y, x) fe = 2, Z - 1 (88) 



i=0 



;-2 



F (^"V, x) J] [1 + F (^y, x)] = y (89) 

i=0 

Clearly, the last equation is the x-deformation of the equivalent equation (28) of 
section 2. We know that, for x close enough to 1, we have I solutions Fj (y,x) 
which have the property that Fj (y, 1) = Vj (y,l) described in section 2 below 
(29). 

From (82-83) we have 

l-l 1-2 

( y ,x) = i+Y, w i (y> x ) = n i 1 + F ( xl y> x )] ( 9 °) 

i=l i=0 

so that 

d+l-2 

*d(y,x)= J] [l + F(x l i/,x)] (91) 

i=0 

The partition function is obtained from the linear system (83) (I functions 
Z d _i (x l y,x) and I solutions of (89) for (y,x) and Wi (y,x)); consequently, 

l-l d+l-2 

Z(x,y,l,l) = J2^(y^) II [^ J rF j (x i y,x)] (92) 

j=0 i=0 

where Aj (y, x) is the ratio of two determinants defined by the linear system (83) 

and where the index j refer to the j th solution Fj (y,x) of equation (89). We do 
not wish to explore further the properties of these determinants (see appendix 
B); let us simply remind to the reader that from (29) in section 2, we know that 
for x — 1, 



We now consider the thermodynamic large <i-limit of the partition function 
with x = — > 1. As before, for x close enough to 1 and y small enough, we 
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may consider the solution F (y, x) to be dominant in the limit. We must use the 
following expansion 



\n[l + F(y,x)]=J2a n (x) y n 



n=l 



together with the property 

" y dv 
/„ v 



f 

Jo 



y_ 

n 



In [l + F(v,x)] = J2 

&n \X) 

n=l 

We may now make a ^ expansion of nf=o~ 2 [l + F £ 3 



and write 



Zth (t,y,l,i) 



1 + V(y,l)\ 3 * [1 + V(£y,l)] 1 



l + V(£y,l)J l + lV(y,l) 
d 



exp [$ (£y) - $ 



. exp 



where 



*(v) 



d f v dv 

/ — \n[l + F (v,x) 
ox J V 



x=l 



(94) 



(95) 



(96) 



(97) 



This ends the calculation of the finite ci-size corrections in the case m — 1. At 
the end of subsection 4.4, we show that the integral in (96) can be expressed in 
terms of the functions Li + , Li 2 or of the Rogers dilogarithm function. 



4.4 The partition function Z(x,y, /,m) 



The partition function in the general case (for any I and m) is calculated in 
appendix B. We simply give the result and its thermodynamic limit. We obtain 



Z r (x,y,l,m) = 



x r d+ dr=ll yr sup(m,0-l 



d+r-2 



m 



*i(v,x) n [i + *;■(*'!/,*)]. 



r-1 



•nN 

i=0 



3=0 



i=Q 



(98) 



In this expression, Aj (y, x) is the ratio of two determinants and is explained in 
Appendix B; the functions Fj (y, x) satisfy the equations 



m— 1 



m—l—1 



x -{m-l) p ^ x ± y ^ x j = x ~ (m ¥ l 1} y ™ | J ^ + F ^ m > z (gg) 

i=0 i=0 
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m—1 



l—m—1 



Yl F (x l ™ +l y,x"j Yl I 1 + F (x™y,x) 

i=0 i=0 



(m-l)fi-l) 

x~- 2 y m I > m (100) 



where Yio 1 [•••] — 1 by convention. The equations (98), (99) and(100) are the x- 
deformed equations (A18), (Al5a) and (Al5b) of appendix A. By identification 
of (98) when x = 1 with (A18), we see that 



Fj (v, i) = v. 



(101) 



A, (y, i) 



m(l + Vj) 
m + lVi 



(102) 



We are now taking the thermodynamic limit when x = £ a — * 1 and when 
the number of states d — > oo. The technique is exactly similar to the one used in 
subsection 4.3; we obtain 



Z r ,th(t,y,l,m) = [l + V(£y,l/m)] 

d r y dv 



h ( 1 + V (y, l/m) \ 3/2 exp [$ (&/) - $ (y)] 



. exp 



lnf 



l + ^(£y,V m )/ m + lV(y,l/m) 
J V ^-ln[l + V (v, l/m)]) + O (103) 



where 



$(y) 



/ -ln[l + FoKx) 
Jo v 



d_ f y dv 

dx 



(104) 



x=l 



Finally, the partition function in the thermodynamic limit is obtained from (103) 
by summing r from 1 to m 



Z th (£,y,l,m) = 



[1 + V (£y, l/m)} 1 - 1 fl + V(y, l/m) \ 3/2 exp [$ (&) - $ (y)} 



1-[1 + V(£y,l /m)] ~ l/m V + V (&/, l/m) J m + lV(y, l/m) 

d T f % ^ [1 + V (v, l/m)]) + O (105) 



. exp 



Using (37), the integral in (105) can be transformed into Li + and Li 2 functions 
since 

/ — \n[l + V(v,l/m)] = Lt + [V{y,l/m)} + — —In 2 [1 + V (y, I /m)ll06) 
Jo v ^ m 



Lin 



V(y,l/m) 
l + V{y,l/m) 



+ ^ In 2 [1 + V (y, l/m)] 

(107) 
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and finally using the equation 



i — JTl 

V(l + V)~ = y 



(108) 



we also have 

r y dv 

Jo 



V 



In [1 + V (v,l/m)) = L 



V (y, l/m) 
l + V(y,l/m) 



+ -lny . \n[l + V(y,l/m)} 



where the function L (y) is the Rogers dilogarithm function 

1 



L(y) = - l 



dv 



\nv In (1 — v) 
+ — 



1-v 



(109) 



(110) 



This result has been obtained by Berkovich and McCoy' 12 ' . The difference between 
their result and ours in (105) is: 

1°) their result is obtained using the steepest descent method and provides 
only (109) while we get the finite d-size corrections, 

2°) they calculated the limit x — > 1 [1 + F(x l y,x)] while we calculated 

the limit d — > oo Yli=o 1 + (C^y,^^ which generates the ^ dependant sub- 
tracted terms. 



In our limit, the average number of particles per state is found to be 

ln ( i+v(y,g) A 

111 \l+V(£y,g)J 



n(t,y) 



inr 1 



111) 



where g = l/m. 

From conformal field theory, we know that the central charge is given by 



6 f 1 dv 



7T 



ln[l + y (y,l/m)] 



l+V(l,g) 



L(l) 



(112) 



Following Hikami' 11 ] we define 



and we obtain the central charge as 



V(l,g) 



l + V(l,g) 



(113) 



H) 



L(l) 



(114) 
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where £ H satisfy 



(115) 
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5 Appendix A 



1°) The function V (y, a) 

We define the function V (y, a) as the solution of the equation 

z{l + z) a ~ l =y ae R (Al) 

which has a small y expansion of the type 

V(y,a)=y+(l-a) y 2 + O (y 3 ) (A2) 

We may calculate the complete y expansion of V k (y, a) for any integer k from 
the contour integral 

,rk, x 1 f , z k (1 + z) a ' 2 (1 + az) 
lm Jc v z{l + z) -y 

where the numerator is z k times the derivative of the denominator and where 
the contour Cy surrounds the point z = V (y, a) . In order to obtain the formal 
y expansion, we deform the contour to include the point z = but exclude the 
point z — — 1. We easily obtain 

V k (y,a) = y k + k (1 - a) y k+l + 

+k(a-i)±(-r r r ^+M^-i)] (116) 

V ' ' r \na + k(a- 1 )-n + l \ n\ V ; 

71=2 J 
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This expansion is in fact valid for any k G R and for instance for k — — 1, we 
obtain the y expansion (5) for the function W(y,g) = V~ l (y,g) which enters 
the expression (3) for the average number of particles per site. From the above 
expansion and from the equation (Al), it is easy to obtain some other useful 
expansions like 

/-, t 7" / xxfe n , , / \ na — k — 1] ... [not — k — n + 1] „ 
1 + V (y, a)f = l + ky-kj_] (-) n l - J— L l - y n 

(A5) 



ni 

n=2 



and 



, / . / , n (na — 1) ... (no — n + 1) „ . . „ s 
In (1 + V (y, a))=y-J2 H " " 1/" (A6) 



71! 

n=2 



We now consider the solutions of the polynomial equation 

z m = (1 + z) m ~ l y m m >l (A7) 

where m and / are non negative integers. The above polynomial in z has m roots 
which are 

Z j = V (q j y,l/m) j = 0, ...,m - 1 (A8) 
where g = exp {2m /m) . 

We now consider the solutions of the polynomial equation 

z m (1 + z) l - m = y m l>m (A9) 

where m and Z are non negative integers. The above polynomial in z has / roots, 
m of which — > when y — > 0; they are 

z . = V (q j y,l/m) j = 0,...,m-l (A10) 

where q = exp (2m /m) . The (/ — m) remaining roots — > —1 when y — > 0. If we 
write z = — 1 — t, we obtain for t a similar equation as (A9). Consequently the 
remaining roots can be written as 

zj = -l-V (q 2 i +1 yr^, -L-) j = m, / - 1 (All) 

where Q = exp (in/ (I — m) ) . 
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2°) Calculation of the partition function Z(y,l,m). 
We define 



Z r (y,l,m) = 5 r , m + CSS^O+r-i yPm+V ( A12 ) 

p=0 

Using the contour integral representation (26) for the combinatoric symbol, we 
write 

\ l J , „ /-, \d+p{m-l)+r-l 

Z r (y, I, m) = S r , m + dz [ - y^ (A13) 

p=0 J °o 



where C surrounds the point z = 0. The sum over p can be extended to oo since 
the contour integral gives for p > E (^) . We obtain 

r r m-r-1 Q , z ) d + r ~ 1 

Z r (yJ,m) = ^-<p dz Z - (A14) 



2wr J c z ™_( 1 + z y 



The convergence of the series requires \z\ m |1 + z\ l ~ m > y which is certainly 
true for large \z\ . Consequently, the contour C^ m surrounds all the roots of the 
polynomial equation 

\m—l 



1 + z) m ~ l y m = m > I (Al5a) 



or 



In (A14), the extra pole at z = when r — m has been taken into account to 
cancel the term 5 r>m in (A13); in that case, the contour C^ m does not surround 
the point z = 0. The equation (Al5a) has m roots 

y. = V (qiy, l/m) j = 0, m - 1 (A16) 

The equation (Al5b) has / roots 

Vj = V(qiy,l/m) j = 0,...,m-l 

Vj = -l-V^Q^y^,-^^ j=m,...,l-l (117) 

with q = exp (2m /m) and Q = exp (in/ (I — m)) . The contour integral (A14) is 
now calculated by summing over all roots 

Sup(rn,l)-1 d+r r 

Z r (y,l,m)= Y K {\ r It (A18) 
j=0 J j 
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Finally the partition function is obtained by summing the geometrical series in 
(A18) from r = 1 to r = m : 

Z(y,l,m)= Y [ -4 — { 3 L r (A19) 

where k (j) is an integer equal to j for j = 0, ...,m — 1 but depends upon the 
definition of (1 + Vj) 1/m in (A17) for I > m and j = m, I - 1. 



6 Appendix B 



In this appendix, we calculate the partition function 

e( 



' < kzl % 



Z r . (x, y, /, m) = Z d , (y, x) = 5 r , m + £ ^^W)' C ^ m _ l)+r _ x (x) y*** 

(Bl) 

for any value of I and m. The method used here is the same as in section 4.3 
for m — 1; we calculate the difference Z d ^ r (y,x) — Z d -\ )T (xy,x) and we use the 
relation (79) between the combinatoric symbols C. The result depends on the 
value of r : 

Z d ,r (y,x) = Z d _ 1>r (xy,x) +x~ {r ~ l) ™ y Z d<T _ x (x^yj r = 2,3, ...,m (B2) 



Z d ,i (y, x) = Z d _ 1A (xy, x)+y Z d ^ m [x™yj (B3) 

In order to use these recursions we may introduce a generalized Fourier transform 
which diagonalizes the system over r and in the same time performs subtractions 
over the number of states d. As a matter of fact, the subtractions over d are 
rational so that we must consider the existence of Z d ^ r (y, x) for any rational d. 
We define 

m 

T djTl (y,x) = VV r x~ LJJ ^ Z d _, r _ 1)± (y,x) n = 0, ...,m-l (B4) 

r=l 

where q = exp(2in / m) . The generalized Fourier transform (B4) can be inverted 

as 

r(m-r)l m —\ 

Z d Ay,x) = ^^y,q- nr T d+{r _ 1)±!n (y,x) r = l,...,m (B5) 

n=0 
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Then, the recursions (B2) and (B3) are equivalent to the following recursion 

T d ,n (y, x) = T d _i, n (xy, x) + q n x" 1 ^ y T d _^ n (x^y, x^j (B6) 

The solution of the above recursion is possible if we introduce two coprime integers 
/' and m! such that — = -^7. It is important all along the solution to keep track 
of the q n terms so that we may at the end perform the inverse Fourier transform 
(B5). Let us decompose n in a unique way (Bezout decomposition) as 

Tfl 

n = n' + sl (mod m) < n' < — -1 < s < m' - 1 (B7) 

m 

Now, the solution of the above recursion depends upon the relative values of / 
and m. 

l°)m>l 

We define 

rn'-l 

V d , n (y, x) = T dtn (y, x)+J2 W i,* {y, ^-^,n (x^V, x) (B8) 

i=i m 

where Q = exp (2m /m!) and where the functions W^ n i (y,x) are determined by 
the condition 

^> dtn (y, x) = Q s U n , (y, x) %-A,, n (x^y, x^j (B9) 

which implies 

m'-l 

^d,n(y,x)= Yl u n' (x^y,xj Vd-i,n(xy,x) (BIO) 

i=0 

The constraints (B8) and (B9) gives the relations 

W ky (y, x) = U n > (y,x) W k ^ n > (x^y,x) k = 2, V - 1, V + 1, m' - 1 
W 1>n ,(y,x) = U n ,(y,x) ifl'^l (118) 

q n ' x~ (J2 ^ 11 y + W V) n' (y, x) = U n > (y, x) W V - X , n ' (x^y, x^j (Bllb) 
where eventually Wo, n ' (y,x) = 1, 

1 = U n > (y,x) W m '-i, n ' (x^y,x) (Bile) 
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These equations can be solved recursively; we obtain 

k-i 



W ky (y, x) = Yl U n , (x^y, xj k = 1, I' - 1 

i=0 



(B12a) 



k-l'-l 



w. 



'k,n'(y,x)= Yl U n' (x™'V,x} W Vin > (x m' y,xj k = I' + 1, m' - 1 

(B12b) 



i=0 



i'-i 



/ (m — l)l "| — |" / j \ 

g n x 2- y + W v , n > (y, x) = [ [ U n > (x^y, xj (B12c) 

i=0 



m'-l'-l 

1 = JJ U n > (x^y,xj Wy >n > (x m ^ L y,x^j 

i=0 



(B12d) 



m — Z 

If we transform (B12c) by changing y — > x ™' y and if we multiply by 
nHo~ Z ~* ^n' ^ x " 7 l/5 ;r j we obtain from (B12d) 



m'-l'-l 



q n ' x (m 2m 1)! y JJ [/ n , (x^y,x) = x ^ F n ,{y,x) (B13) 

i=0 

where 1 + F n /(y,x) = YYILq 1 {^^y-i^ ■ Again, if we multiply the above 
equation by (mf — 1) similar equations, each being deduced from the precedent 
one by the transformation y — > x^y, we obtain 



m'-l 

° n f » 

i=0 



/ I X 



g x 



, , [(m-l)l'-m'+l] 

,n m „— — — 



m'-l'-l 



y m ' J] [l + F n , (x^y,x)] (119) 

i=0 



If we multiply (B14a) by m/m' different equations, each being deduced from the 
precedent one by the transformation y — > x™y, we obtain 



m—l 



X 



-(m—l) 



i=0 



(m-l)(i-l) 
X 2 



m—l—1 



y m J] [l + F n , (x^y,x)] (B14b) 

i=0 
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which is the x-deformed equation (Al5a). At x — 1, equation (B14a) is polyno- 
mial of degree to'. For a given n', we have to' solutions F^,\ j = 1, ...,to'; this 
analysis remain true by x deformation for x not too far from 1. On the other 
hand, equation (B14b) is n' independant and has to solutions F u (y, x) which 
are the same globally as the solutions F^p for all n' and all j. Considering for a 
given n', the to' solutions F^,\ j = 1, ...,m', the system (B8) is a linear system 
of m' equations for m' unknown functions T d _j^, (x~^y,x^ , i = 0,...,m' — 1. 
Consequently, we may write 

rn' 

T d _^(x^y,x)=Q~ is J2 A tu)(V' x ) i = 0,...,m'-l (B15) 

i=i 

where A™'^ (y, x) is the ratio of two determinants containing only the functions 

wfy (y, x) , the powers of Q being extracted separately. The functions A"^ (y, x) 
satisfy the relations 

A&) (y,x) = W#>_ 1>B , (x*y,x) ^L W) (x*y,x) i = 1, ...,m' - 1 (B16) 

and the orthogonality relations 

m'-l 

J2w$(y,x) A? {j) (y,x) = 5 



3 



i=Q 
m' 



J2 w tt>(y^) A iij)(y^) = (120) 

Finally, using 

d-l 

(y, x) = n [ x + *)] *) ( B18 ) 

i=0 

the solution of the recursion (B6) is 

m' d— 1 

7^ (x*y,x) = Q~ is J2 A t(J) M II [ l + ^ (* d V,*) 
™ i=i fc=o 

(B19) 

We are now in position to calculate the inverse Fourier transform (B5). We write 
(r — 1)/' = Km! — h with K integer and < h < to'; then, we find 

r(m — r)l 

n',j 

d+K-1 

. jj [l+^p-^^j;'^,.) (m) 
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where the sums over n' and j runs respectively from to ^ — 1 and from 1 to 
w!\ in (B20), we introduced the functions 



m'—l 

*$(y,x) (B2i) 

s=0 

which may be determined from the values of Z d r (y, x) in (Bl) for < d < I + 1. 
We find 

$ { i)(y,x) = W l 9l(y,x) (B22) 

Now, the expression (B20) can be transformed using successively (B16), (Bile), 
(B22), (B12d) and the expression of [1 + F n i] in terms of U n >; we obtain 

m r(m— r)l d+r—2 



Z*,Av,x) = ^f^E^Xo-)^^) II h + F^(x k y,x)}.(122) 

lit L 

n ,j k=0 
r(m'-l')-l 

■ n K'^H 



i=0 

Finally, using (B13) r times and (Bile), we get 

rd I r ( r ~ 1 ) r m— 1 d+r-2 

Z d , r (y,x) = X " y E A "^' X ) II [1 + K(x k y,x)]. 

u=0 k=0 

■n[^«(^" 1+(r " i) "y^)]" 1 ( 123 ) 

i=0 

where u = n' + (j — 1) Equation (B24) is the x deformation of (A18) for 
m > I. 

2°)l>m: 

We define 

V-i 

tt d , n (y, x) = T d , n (y, x) + J2 Q* x ) T d-±r,n (x^y, x) (B25) 

i=i m 

where Q = exp(2in / W) . The functions Wi j7l / (y,x) are constrained by the condi- 
tions 

* d>n (y, x) = Q s WV (y, x) * d _^ in (x^y, x) (B26) 

which implies 
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m'-l 



^d,n(y,x) = Y[ W iy (x^y,x^ V d _ hn (xy,x) (B27) 

i=0 



We obtain the following relations 

W k ,n'(y,x) = W 1>n > (y,x) W k -i, n ' (x^y,xj 
fork = 2,...,m'-l,m , + l,...,Z'-l 

1 + W m/ , a/ (y, x) = W 1>n , (y, x) W m >_ hn , (x^y, xj 

q n ' x"^ 1 - y = W hn ,(y,x) W V -i, n > (x^y,x) 

These three equations are easily transformed into 

fc-i 

W k , n ' (y, x ) = Y\. Wl > n ' ( x ^> x ) k = l, m'-l 



(124) 
(B28b) 

(B28c) 
(B29a) 



i=0 



fc— m' — 1 



W k , n >(y,x) = Yl W hn' (x^y,xj W m ,, a , (x^y^x^j k = rri + 1, 1 

(B29b) 



i=0 



m'-l 



1 + 



W m ', n ' (y,x) = Yl w i,n' (x^y,xj 

i=0 



(B29c) 



, (m-l)i 

q x y 



/'-m'-l 



n Wi,„' W m / jn / {xh^y,x^ (B29d) 

i=0 



If we multiply this last equation by (m' — 1) other equations, each being trans- 
formed from the one before by making y — > x^y, we obtain 



m'-l 



i'-m'-l 



Yl Fw [xH^y.x^ J] [l + F n , (x^y,x)\ = q n ' m ' x 

i=0 i=0 



(m-l)i' 



(B30a) 



where F n i(y,x) = W m > jn r (y,x) . Again, if we multiply equation (B30a) by m/m' 
different equations, each being transformed from the one before by making y — > 
x™y, we obtain 



m— 1 



I— m— 1 



H F n , (x^y, x) H [l + F n , (x*y, x)] = aT 1 " 11 "^ 11 y m (B30b) 

i=0 i=0 
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which is the x-deformed equation (Al5b). At x — 1, equation (B30a) is polyno- 
mial in F n r of degree I'. For a given n', we have /' solutions; this analysis remains 
true by x deformation and for x not too far from 1. Considering the /' solutions 
Fffi, j = the system (B25) is a linear system of /' equations for I' 

unknown functions T d ^ (^x^y, x^j , i — 0, V — 1. We may write 

V 

T d -^, n (x^y,x) =Q- isJ £ A t(j)(y> x ) < = 0,...,Z / -1 (B31) 

where A™',^ (y, x) is the ratio of two determinants containing only W- ^, (y, x) 
functions, the powers of Q being extracted separetely. It is easy to convince 
oneself that 

(x*V> *) = W i% (I/. *) ^ x ) ^ = 0, 2 (B32) 

and that we have the following orthogonality relations 



i=0 

v 

52w$(y,x) At U) (y,x) = 5 hk (125) 
j'=i 



Finally, using (B27) and (B29c), we have 

d-i 

*2i & x ) = II [ X + 0^> *)] *) (B34) 

i=0 

so that the solution of the recursion (B6) is 

«' d-l 

i=i fc=o 

(B35) 

for i = 0, -1. 

We must now perform the inverse generalized Fourier transform (B5). Again, 
we write (r — 1)1' = Km' — h where K is an integer and < h < m'. Then 
equation (B20), where j is summed from 1 to I', and equation (B21) are valid. 
Again, the initial conditions that is the values of Zd, r (y, x) from (Bl) for < d < 

I + 1 determine the functions 

l'-m'-l 

*$(V,X)= II W$(x*y,x) (B36) 

i=0 
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Using now (B29c), the definition of F$ (y,x), (B32) and (B36), we may write 



Zd,r (y, x) 



X 



r(m—r)l 
2m 

m 



d+r-2 



II [^ + F^(x k y,x)_ 

n',j k=0 



r(/'-m')-l 



n < 



C?) („d+r-l+- 



i=0 



Finally, using (B29d) r times, we get 



z d) r (y, x) = 



X 



rd+ 



d+r-2 



(126) 



m 



J2 A u(y,x) fl [1 + F u (x k y,x)]. 



r-1 



u=0 



k=0 



.H^fx^-^-^x)]' 1 (127) 

i=0 

where u = n' + (j — 1) y. Equation (B38) is the x deformation of (A18) for I > m. 
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Forbidden configuration for m=3 
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Fig. 1: "Composite particles' 
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